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Abstract 

For a successful detection of gravitational waves by LISA, it is essential to construct theoretical waveforms in a 
reliable manner. We discuss gravitational waves from an extreme mass ratio binary system which is expected to be a 
promising target of the LISA project. 

The extreme mass ratio binary is a binary system of a supermassive black hole and a stellar mass compact object. 
As the supermassive black hole dominates the gravitational field of the system, we suppose that the system might 
be well approximated by a metric perturbation of a Kerr black hole. We discuss a recent theoretical progress in 
calculating the waveforms from such a system. 

I. INTRODUCTION 

We consider a method to construct theoretical templates of gravitational waves LISA will detect. We are planning 
to operate LISA around a few years, and we may observe gravitational waves for this time scale. One of promising 
targets for LISA is a binary system of a supermassive black hole and an inspiralling compact object, and we discuss a 
method to calculate a gravitational waveform coming from such a system. Since the supermassive black hole dominates 
the gravitational field of the system, we expect that a linear metric perturbation of a Kerr black hole induced by a 
point particle might be a useful tool for this calculation. 

Because we already have a method to calculate a metric perturbation of a given orbit of the particle, it was widely 
believed that the key issue here is to calculate a self-force, so that one can derive the orbit of the particle. We derived 
a formal expression of the self-force [1] and we come to have a promising method to calculate the self-force to the 
particle [2]. However, there was a concern on this approach because of a gauge dependence of the self-force. 

We discuss that the self- force is totally gauge dependent. By a special choice of the gauge condition, the self- force 
vanishes for the entire time interval where the linear metric perturbation is valid. (See Sec. II.) On the other hand, 
gravitational radiation reaction is observable by asymptotic gravitational waves, and we expect that a physically 
reasonable self-force must describe this effect. Thus, the vanishing self-force suggests that we need a substantial 
extension of the theoretical basis in calculating the self-force, especially, of the gauge condition. 

We find a key idea to solve this problem is an adiabatic approximation of the metric perturbation, which is 
different from a usual metric perturbation scheme. The idea of an adiabatic approximation of the metric changes our 
understanding of the orbital evolution and the validity of the self-force was clarified for the first time. We find that 
the class of gauge conditions used in the radiation reaction formalism in Ref. [3] is compatible with the adiabatic 
approximation of the metric perturbation. With a further choice of the gauge condition, we find that the radiation 
reaction formula might describe the orbit sufficiently long enough for the LISA project. 

II. SELF-FORCE? A MYTH? 

Gravitational radiation reaction is physically real since we can define the gravitational flux at the asymptotically 
flat region, and we expect to see this effect on the orbit by a modulation of an observed gravitational waveform. In 
order to predict this modulation, we were motivated to calculate the self-force for the belief that the self-force tells us 
the gravitational radiation reaction effect on the orbit. However, we find that the self-force could vanish consistently 
with a usual perturbation method irrespective to gravitational radiation reaction. 
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For an approximation of an extreme mass ratio binary system, we take the mass ratio of the binary /i as a small 
parameter for the expansion. In a usual perturbation scheme, we expand the metric and the orbit from background 
quantities of the supermassive black hole as 

9* P = giT + M$ + + ■ ■ ■ . (2-1) 

z a {\) = z ^ kerr)a {\) + iiz {1)a {\) + M 2 z (2) "(A) + • • • , (2.2) 

where g^ r is a Kerr metric, and z( feerr ) a (A) is a geodesic around it with A as an orbital parameter. For a valid 
perturbation, we assume 

0(1) -9% >M$ > A$ > (2.3) 

0(1) ~ z^ a (X) > M zW«(A) > M 2 z( 2 )"(A) > ••• , (2.4) 

and the perturbative equations are derived by expanding the Einstein equation with respect to /i. 

One can see that this perturbation scheme allows only a small deviation of the orbit from the background geodesic. 
Because the orbit will deviate from the background geodesic eventually by gravitational radiation reaction, this 
perturbation scheme is valid in a finite time interval. 

We consider a gauge transformation of the orbit in this finite time interval where the perturbation is valid. Because 
the deviation from a geodesic is 0(/i), one can always take a gauge transformation to eliminate the deviation during 
this whole time interval, i.e. = 0, thus, the self-force totally vanishes consistently with the usual 

perturbation scheme. 



III. ADIABATIC METRIC PERTURBATION 



The reason for this gauge problem of the self-force is that the perturbation is valid only in a short time scale, so 
that we can only describe a small deviation of the orbit from the background geodesic. In order to solve this problem, 
we propose an extension of the perturbation scheme. 

Bound geodesies around a Kerr black hole with a convenient orbital parameter A are characterized by three constants 
of motion £ a , a = E,L,C, and four phase constants A fc , b = r 7 9, C c , c — t, (j) as in Ref. [3], and we denote them 
collectively by 7 = {£ a , X b , C c }. From the result of Ref. [3], the orbital deviation becomes 0(1) when A ~ 0(^i _1 / 2 ). 
This time scale is called the dephasing time. From the discussion in Sec. II, the usual perturbation scheme is valid 
only in this time scale. Our purpose here is to extend this time scale of validity by modifying the metric perturbation 
scheme. 

In the usual perturbation scheme, the linear metric perturbation is induced by a geodesic, and we can write it as a 
function of 7, h a p{x) = h a p{^; x). Here we use a class of gauge conditions where the metric perturbation is derived 
by the Green function in the Boyer-Lindquist coordinates as 

h af3 (x) = f dx'G a/3a , (x,x')T a 'P'(x'), (3.1) 

G aSi a , p , (x, x') = G aP a , p , (t-t',<f>-<j>',r,r',6, 0') . (3.2) 

We assume an adiabatic approximation to the orbit as in Ref. [3], and consider the evolution of "geodesic constants" by 
the effect of gravitational radiation reaction. We write the "geodesic constants" as functions of the orbital parameter, 
7(A). We consider to extend this idea to the metric perturbation. We foliate the spacetime by spacelike hypersurfaces. 
Using the fact that the foliation surfaces intersect with the orbit, we define the foliation function f(x) by the orbital 
parameter at the intersection of the orbit and the surfaces as f(z(X)) — A. We define the adiabatic metric perturbation 
on the foliation surface f(x) by the linear metric perturbation induced by the geodesic j(f) as 

h a c $(x) = h af3 (j(f);x). (3.3) 

In order to see the validity of the adiabatic metric perturbation, we operate the linearized Einstein operator and 
we have 

Gi> ad ]=T^[ 7 (/)]+Ai> ad ], (3.4) 

where an extra term A^J appears because the adiabatic metric perturbation is not induced by a geodesic. It is notable 
that T a p[y(f)] is the stress-energy tensor of a point particle moving along the orbit with the effect of gravitational 
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radiation reaction. The adiabatic metric perturbation solves the Einstein equation to the accuracy of 0{fi) as long 
as 0(jjl) > A$ holds. Since A^g depends on the A derivatives of the orbital "constants" 7(A), the validity of the 
adiabatic metric perturbation depends on how the orbit evolves. We find that, under the gauge condition (3.2), the 
adiabatic metric perturbation is valid in the radiation reaction time scale 0(^. _1 ) > A. 

Because the time scale of validity is much longer than that of the linear perturbation, the orbital deviation could 
be O(l) in this metric perturbation scheme and one cannot eliminate the orbital deviation from the geodesic by a 
gauge transformation, thus, the self-force cannot vanish by a gauge transformation and is physically meaningful to 
predict the orbit in this class of gauge conditions. 

IV. ADIABATIC EVOLUTION OF THE ORBIT 

Ref. [3] finds that the infinite time averaged self-force acting on the "constants of motion", < (d/d\)£ a >, can 
be derived by the radiative metric perturbation in a gauge invariant way by a simple symmetry property of a Kerr 
spacetime. It also find that < (d/d\)£ a > makes a dominant contribution to the orbital evolution by a perturbative 
analysis of the orbit. 

The adiabatic metric perturbation is valid in the radiation reaction time scale and, in this time scale, the orbital 
evolution is non perturbative. As for the "constants of motion" , one can still deal with the evolution by a perturbation, 
and we find that < {d/d\)£ a > dominantly determines the evolution of them. The subdominant part of the "constants 
of motion" is O(p). The evolution of the "phase constants" becomes non perturbative beyond the dephasing time, 
however, we find the perturbative results in Ref. [3] are qualitatively correct. The dominant part of their evolution is 
described only by < (d/d\)£ a >. and the subdominant part behaves as O(jiX). 

Because the subdominant part of the evolution of "constants of motion" is smaller than 0(1) in the radiation 
reaction time scale, one can find a gauge condition which can totally eliminate this part. We call this gauge condition 
by the radiation reaction gauge, and, in this gauge, the self-force can be written only by the infinite time averaged 
self-force as 

(d/d\)£ a =< (d/d\)£ a > . (4.1) 

The self-force is expected to have the conservative part as well as the radiative part, however, in this gauge 
condition, the conservative part is integrated out and appears to be the small shift of initial values. We also note 
that the radiation reaction gauge condition applies to an orbit of a spinning particle. It is known that the orbit of 
the spinning test particle slightly deviates from a geodesic [5] because of the coupling of its spin and the background 
curvature. This effect of the deviation can also be absorbed by a small change of the initial values by the radiation 
reaction gauge condition. 

V. SELF-FORCE CALCULATION VS RADIATION REACTION FORMULA 

Here we try to compare these two approaches to calculate the orbit; the regularization calculation of the self-force 
[2] and the radiation reaction formula discussed in Sec. IV. There could be two issues for comparison; 

1) Theoretical issue; whether the resulting orbit is physically acceptable or not 

2) Practical issue; whether the approach is practically available for generating templates for the LISA project 

Gravitational radiation reaction is a physically real effect. Because the momentum flux carried by gravitational 
waves is well defined, it is reasonable to expect that the orbital evolution must be consistent to this effect. There 
is a wide-spread belief that the self-force includes the effect of gravitational radiation reaction, however, this is not 
trivial in the usual perturbation scheme because of the gauge freedom as we argue in Sec. II. We consider that the 
regularization calculation of the self-force [2] is not guaranteed to include this radiation reaction effect in general. On 
the other hand, the radiation reaction formula uses a certain class of gauge conditions as we argue in Sec. Ill, and the 
orbital evolution is described consistently with asymptotic gravitational radiation. 

For a practical issue, an idea is proposed to calculate the metric perturbation in a Kerr black hole available for the 
self-force calculation [6], however, it is not coded yet, and a more crucial problem is whether the resulting self-force 
is consistent with gravitational radiation reaction. Even if it is possible, the numerical convergence of the self-force 
is not optimal since it is a point splitting regularization. On the other hand, the radiation reaction formula is, in the 
end, just an extension of the balance formula, and various numerical codes of the balance formula were already made 
successfully [7] . There is also known a semi-analytic technique which may substantially increase the efficiency of the 
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code [8] . In this approach, we have the most optimal convergence because the regularization is done in a quite natural 
manner, and we consider it helpful to have an efficient calculation. 
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